This code is Script 4d for Kitchel et al. “Biotic homogenization, the exception and not the rule for marine fish communities” manuscript.

SESSION INFO TO DO

library(data.table)
library(vegan)
library(sf)
library(concaveman) #polygon around points
library(betapart) #allows us to partition beta diversity
library(geosphere)
library(ggpubr) #stat_regline_equation
library(nlme)
library(mgcv) #to make gam
library(cowplot)
library(lme4)

#Pull Dissimilarity Means
distances_dissimilarities_allyears.r <- readRDS(here::here("output","dissimilarities", "distances_dissimilarities_allyears.r.rds"))

#make survey and survey unit factors
distances_dissimilarities_allyears.r[,survey:=factor(survey)][,survey_unit:=factor(survey_unit)]

#adjust years
distances_dissimilarities_allyears.r[,year_adj := year-min(year)+1]

#add new variable for year in sequence per region
distances_dissimilarities_allyears.r[,first_year := min(year),.(survey_unit)]
distances_dissimilarities_allyears.r[,last_year := max(year),.(survey_unit)]

#distances_dissimilarities_allyears.r[,year_in_seq := year-first_year+1]

distances_dissimilarities_allyears.r[,years_sampled := last_year-first_year+1]

###Palette for Plotting Palette for plotting all 37 survey units

survey_unit.list <- levels(distances_dissimilarities_allyears.r[,survey_unit])

palette_42 <- c(
  "#5A5156", #AI
  "#DF00DB", #BITS-1
  "#DB8EDA", #BITS-4
  "#F6222E", #CHL
  "#F8A19F", #DFO-NF
  "#16FF32", #DFO-QCS
  "#325A9B", #EBS
  "#3283FE", #EVHOE
  "#FEAF16", #FR-CGFS
  "#fccb6d", #GMEX-Fall
  "#1C8356", #GMEX-Summer
  "#C4451C", #GOA
  "#85660D", #GRL-DE
  "#B0009F", #GSL-N
  "#BF79B8", #GSL-S
  "#1CBE4F", #ICE-GFS
  "#782AB6", #IE-IGFS
  "#90AD1C", #MEDITS
  "#6B003A", #NAM
  "#A75B00", #NEUS-Fall
  "#E3B072", #NEUS-Spring
  "#02E8B6", #NIGFS-1
  "#97E7D5", #NIGFS-4
  "#B00068", #Nor-BTS-3
  "#00B9E3", #NS-IBTS-1
  "#95E2F4", #NS-IBTS-3
  "#B3CE73", #NZ-CHAT
  "#689500", #NZ-ECSI
  "#364d02",#NZ-SUBA
  "#AAF400", #NZ-WCSI
  "#AA0DFE", #PT-IBTS
  "#7f9eb8", #ROCKALL
  "#FA0087", #S-GEORG
  "#DEA0FD", #SCS-Summer
  "#FCEF88", #SEUS-fall
  "#A59405", #SEUS-spring
  "#FCE100", #SEUS-summer
  "#544563", #SWC-IBTS-1
  "#a37fc7", #SWC-IBTS-4
  "#C075A6", #WCANN
  "#BDCDFF", #ZAF-ATL
  "#003EFF"  #ZAF-IND
)

color_link <- data.table(survey_unit = survey_unit.list,hex = palette_42)

Add names for plotting


name_helper <- data.table(Survey_Name_Season = c("Aleutian Islands",
                                    "Baltic Sea Q1",
                                    "Baltic Sea Q4",
                                    "Chile",
                                    "Newfoundland",
                                    "Queen Charlotte Sound",
                                    "Eastern Bering Sea",
                                    "Bay of Biscay",
                                    "English Channel",
                                    "Gulf of Mexico Summer",
                                    "Gulf of Alaska",
                                    "Greenland",
                                    "N Gulf of St. Lawrence",
                                    "S Gulf of St. Lawrence",
                                    "Iceland",
                                    "Irish Sea",
                                    "Mediterranean",
                                    "Namibia",
                                    "NE US Fall",
                                    "NE US Spring",
                                    "N Ireland Q1",
                                    "N Ireland Q4",
                                    "Barents Sea Norway Q3",
                                    "N Sea Q1",
                                    "N Sea Q3",
                                    "Chatham Rise NZ",
                                    "E Coast S Island NZ",
                                    "W Coast S Island NZ",
                                    "Portugal",
                                    "S Georgia",
                                  "Scotian Shelf Summer",
                                  "SE US Fall",
                                  "SE US Spring",
                                  "SE US Summer",
                                  "W Coast US",
                                  "Atlantic Ocean ZA",
                                  "Indian Ocean ZA",
                                   "Rockall Plateau",
                                  "Scotland Shelf Sea Q1",
                                  "Scotland Shelf Sea Q4",
                                  "Falkland Islands",
                                  "Gulf of Mexico Fall",
                                  "Barents Sea Norway Q1",
                                  "Sub-Arctic NZ",
                                  "Scotian Shelf Spring"),
                          survey_unit = c(
                                  "AI",        
                                  "BITS-1",    
                                  "BITS-4",    
                                  "CHL",       
                                  "DFO-NF",    
                                  "DFO-QCS",   
                                  "EBS",       
                                  "EVHOE",     
                                  "FR-CGFS",   
                                  "GMEX-Summer",
                                  "GOA",       
                                  "GRL-DE",    
                                  "GSL-N",     
                                  "GSL-S",     
                                  "ICE-GFS",   
                                  "IE-IGFS",   
                                  "MEDITS",    
                                  "NAM",       
                                  "NEUS-Fall", 
                                  "NEUS-Spring",
                                  "NIGFS-1",   
                                  "NIGFS-4",   
                                  "Nor-BTS-3", 
                                  "NS-IBTS-1", 
                                  "NS-IBTS-3", 
                                  "NZ-CHAT",   
                                  "NZ-ECSI",   
                                  "NZ-WCSI",   
                                  "PT-IBTS",   
                                  "S-GEORG",   
                                  "SCS-SUMMER",
                                  "SEUS-fall", 
                                  "SEUS-spring",
                                  "SEUS-summer",
                                  "WCANN",     
                                  "ZAF-ATL",   
                                  "ZAF-IND",
                                  "ROCKALL",
                                  "SWC-IBTS-1",
                                  "SWC-IBTS-4",
                                  "FALK",
                                  "GMEX-Fall",
                                  "Nor-BTS-1",
                                  "NZ-SUBA",
                                  "SCS-SPRING"
                          ))


color_link <- color_link[name_helper, on = "survey_unit"]

##Make GAMs

Jaccard

bray_curtis_total_gam <- gam(bray_curtis_dissimilarity_total_mean ~ year + s(survey_unit, year, bs = "fs", m = 1),#random smooth
                            data = distances_dissimilarities_allyears.r)

##Make LMERS

Bray These all converged

#running with lme instead of lmer gave same results, but allowed for calculation of p-value
bray_curtis_total_lme <- lme(bray_curtis_dissimilarity_total_mean ~ year_adj, random = (~1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)

#but also run with lmer for confint
bray_curtis_total_lmer <- lmer(bray_curtis_dissimilarity_total_mean ~ year_adj + (1 + year_adj|survey_unit),data = distances_dissimilarities_allyears.r)
Warning: Model failed to converge with max|grad| = 0.137452 (tol = 0.002, component 1)
summary(bray_curtis_total_lme)
Linear mixed-effects model fit by REML
  Data: distances_dissimilarities_allyears.r 

Random effects:
 Formula: ~1 + year_adj | survey_unit
 Structure: General positive-definite, Log-Cholesky parametrization
            StdDev      Corr  
(Intercept) 0.108413092 (Intr)
year_adj    0.001972694 -0.763
Residual    0.023981944       

Fixed effects:  bray_curtis_dissimilarity_total_mean ~ year_adj 
 Correlation: 
         (Intr)
year_adj -0.791

Standardized Within-Group Residuals:
       Min         Q1        Med         Q3        Max 
-4.8791736 -0.4684941  0.0317567  0.5158154  4.3793155 

Number of Observations: 897
Number of Groups: 42 
anova(bray_curtis_total_lme)

bray_curtis_total_coefs <- data.table(coef(bray_curtis_total_lme))
bray_curtis_total_coefs[,survey_unit := rownames(coef(bray_curtis_total_lme))][,Year := round(year_adj,5)][,Intercept := round(`(Intercept)`,2)]
View(bray_curtis_total_coefs)

bray_curtis_total_coefs <- bray_curtis_total_coefs[color_link, on = "survey_unit"]

bray_curtis_total_coefs.exp <- bray_curtis_total_coefs[,.(Survey_Name_Season, Intercept, Year)]

#export this table
fwrite(bray_curtis_total_coefs.exp, file = here::here("output","bray_curtis_total_coefs.exp.csv"))

Get LMER model as predictions


# need to sort out year in seq versus overall year models
#new data for lmer
lmer_year <- seq(min(distances_dissimilarities_allyears.r[,year]), max(distances_dissimilarities_allyears.r[,year]), by = 0.1)

lmer_year_adj <- seq(min(distances_dissimilarities_allyears.r[,year_adj]), max(distances_dissimilarities_allyears.r[,year_adj]), by = 0.1)

#predict average lmer
lmer_bray_total_predictions <- data.table(year = lmer_year, year_adj = lmer_year_adj)

#confidence intervals
bray_curtis_total_lmer_confint <- confint(bray_curtis_total_lmer)
Computing profile confidence intervals ...
#populate data table of lmer predictions
lmer_bray_total_predictions[,bray_curtis_lmer_preds := fixef(bray_curtis_total_lmer)[[1]] + fixef(bray_curtis_total_lmer)[[2]] * year_adj][,bray_curtis_lmer_preds_lowerCI := bray_curtis_total_lmer_confint[5] + bray_curtis_total_lmer_confint[6] * year_adj][,bray_curtis_lmer_preds_upperCI := bray_curtis_total_lmer_confint[11] + bray_curtis_total_lmer_confint[12] * year_adj]

###Move forward with Jaccard Total (others to supplement)

Coefficients for LMER by survey_unit

#unique survey unit year combos
survey_unit_sampling_years <- unique(distances_dissimilarities_allyears.r[,.(survey_unit, year_adj, year, years_sampled)])

# see group coefficients
model_coefs_reduced <- data.table(transform(as.data.frame(ranef(bray_curtis_total_lmer)), lwr = condval - 1.96*condsd, upr = condval + 1.96*condsd))
#https://stackoverflow.com/questions/69805532/extract-the-confidence-intervals-of-lmer-random-effects-plotted-with-dotplotra


#ONLY SLOPES
model_coefs_reduced <- model_coefs_reduced[term == "year_adj",]

model_coefs_reduced[,survey_unit := grp][,year_adj := condval]

#merge with duration of survey
model_coefs_reduced_length <- model_coefs_reduced[survey_unit_sampling_years, on = "survey_unit"]


model_coefs_reduced_length[,years_sampled := as.numeric(years_sampled)][,Directional_Change := ifelse(year_adj > 0, "Differentiation","Homogenization")]

#does it cross zero?
model_coefs_reduced_length[,significant := ifelse(lwr >0 & upr>0,T,ifelse(lwr<0 & upr<0,T,F))]

#significant directional change
model_coefs_reduced_length[,Directional_Change_sig := ifelse(year_adj > 0 & significant == T, "Differentiation",ifelse(year_adj < 0 & significant == T, "Homogenization", "No trend in\ndissimilarity"))]


#min max distances_dissimilarities
min_bray_reduced <- min(distances_dissimilarities_allyears.r[,bray_curtis_dissimilarity_total_mean], na.rm = T)
max_bray_reduced <- max(distances_dissimilarities_allyears.r[,bray_curtis_dissimilarity_total_mean], na.rm = T)

model_coefs_reduced_length <- model_coefs_reduced_length[color_link, on = "survey_unit"]

#delete any NAs
model_coefs_reduced_length <- na.omit(model_coefs_reduced_length, cols = "significant")

#order table by coefficient
setorder(model_coefs_reduced_length, year_adj)

jaccard_total_model_coefs_reduced_length.unique <- unique(model_coefs_reduced_length[,.(condval,condsd, lwr, upr, survey_unit, year_adj, years_sampled, Directional_Change, hex, Survey_Name_Season, significant, Directional_Change_sig)]) 

#extract color hexes
#year adj coef order
color_year_adj_order <- jaccard_total_model_coefs_reduced_length.unique[,hex]

#alphabetical order
jaccard_total_model_coefs_reduced_length.unique.alpha <- setorder(jaccard_total_model_coefs_reduced_length.unique, Survey_Name_Season)

jaccard_total_model_coefs_reduced_length.unique.alpha[,trend_color := ifelse(Directional_Change_sig == "Homogenization", "#e7ac5b", ifelse(Directional_Change_sig == "Differentiation","#91c874","#cbbfde"))]

color_alpha_order <- jaccard_total_model_coefs_reduced_length.unique.alpha[,hex]
color_alpha_order_bytrend <- jaccard_total_model_coefs_reduced_length.unique.alpha[, trend_color]

saveRDS(jaccard_total_model_coefs_reduced_length.unique, here::here("output","region_stats","jaccard_total_model_coefs_reduced_length.unique.Rds"))

Bar Plot Coefficient LMER

jaccard_TOTAL_Dissimilarity_Coef_errorbar_reduced <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, fill = Directional_Change_sig, color = Directional_Change_sig), stat = 'identity', shape = 21) +
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #total jaccard dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(colour = color_year_adj_order, face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.2,0.8), legend.direction = "vertical", legend.background = element_rect(fill = "transparent"))
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: labelWarning: Vectorized input to `element_text()` is not officially supported.
ℹ Results may be unexpected or may change in future versions of ggplot2.
#pull legend for homogenization
directional_change_legend_plot <- jaccard_TOTAL_Dissimilarity_Coef_errorbar_reduced + 
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend") +
  scale_size_binned(range = c(1,8), name = "Survey period length", guide = "none") +
  theme(legend.position = "right", legend.background = element_rect(fill= "transparent")) +
  guides(colour = guide_legend(override.aes = list(size=6)))
Scale for fill is already present.
Adding another scale for fill, which will replace the existing scale.Scale for colour is already present.
Adding another scale for colour, which will replace the existing scale.Scale for size is already present.
Adding another scale for size, which will replace the existing scale.
jaccard_TOTAL_Dissimilarity_Coef_errorbar_reduced


#ALT grey scale
jaccard_TOTAL_Dissimilarity_Coef_errorbar_reduced_greyscale <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, fill = Directional_Change_sig, color = Directional_Change_sig), stat = 'identity', shape = 21) +
  scale_fill_manual(values = c("white","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_color_manual(values = c("black","black","grey"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #total jaccard dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.25,0.8), legend.direction = "vertical")
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: label

Alternatively, we color this plot by trend experienced

jaccard_total_Dissimilarity_Coef_errorbar_reduced_colorbytrend <- ggplot() +
    geom_errorbar(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, ymin = lwr, ymax = upr), fill = "grey", width = 0) + #add confidence intervals
  geom_point(data = model_coefs_reduced_length, aes(x = reorder(Survey_Name_Season, year_adj) , y = year_adj, label = Survey_Name_Season, size = years_sampled, color = Directional_Change_sig), stat = 'identity') +
  scale_color_manual(values = c("#73BA4D","#E0962C","#cbbfde"), name = "Dissimilarity trend", guide="none") +
  scale_size_binned(range = c(1,8), name = "Survey period length\n(years)") +
  geom_hline(yintercept = 0) +
  scale_y_continuous(breaks = seq(-0.005, 0.0075, by = 0.0025), labels = c("-0.005","","0", "", "0.005",  "")) +
  ylim(-0.015,0.006) +
  xlab("Survey unit") +
  ylab("β-diversity trend") + #total jaccard dissimilarity trend
  coord_flip() +
  theme_classic() +
  theme(axis.text.y = element_text(face = "bold"), axis.title.y = element_blank(), axis.text.x = element_text(size = 15), axis.title.x = element_text(size = 15), legend.position = c(0.25,0.8), legend.direction = "vertical", legend.text = element_text(size = 14), legend.title = element_text(size = 15), legend.background = element_rect(fill = "transparent"))
Warning: Ignoring unknown parameters: `fill`Warning: Ignoring unknown aesthetics: labelWarning: Ignoring unknown aesthetics: labelScale for y is already present.
Adding another scale for y, which will replace the existing scale.

Wavy Line Plot for GAMs

Generate predicted values


#add colors and names to full dissimilarity data table
distances_dissimilarities_allyears.r <- distances_dissimilarities_allyears.r[color_link, on = "survey_unit"]

#generate new data to smooth lines (need year and season survey combinations)
year_survey_unit_expand.dt <- data.table(survey_unit = as.character(NULL), year = as.numeric(NULL), year_adj = as.numeric(NULL ))

for (i in 1:length(survey_unit.list)) {
  #generate year vectors
  survey_unit_years <- unique(distances_dissimilarities_allyears.r[survey_unit == survey_unit.list[i],.(survey_unit, year, year_adj)])
  
  years <- seq(min(survey_unit_years[,year]), max(survey_unit_years[,year]), by = 0.1)
  
  year_adjust <- seq(min(survey_unit_years[,year_adj]), max(survey_unit_years[,year_adj]), by = 0.1)
  
  year_survey_unit_expand.dt_addition <- data.table(survey_unit = survey_unit.list[i], year = years, year_adj = year_adjust)
  
  year_survey_unit_expand.dt <- rbind(year_survey_unit_expand.dt, year_survey_unit_expand.dt_addition)
}

#add colors and names to full year and survey unit combination table
year_survey_unit_expand.dt <- year_survey_unit_expand.dt[color_link, on = "survey_unit"]

Get model as predictions

#for plotting, get model as predictions
bray_curtis_total_gam_predictions <- predict(bray_curtis_total_gam, se.fit = TRUE, newdata = year_survey_unit_expand.dt)

#merge into table
year_survey_unit_expand.dt[,bray_glm_mod_fit := bray_curtis_total_gam_predictions$fit][,bray_glm_mod_fit_SE := bray_curtis_total_gam_predictions$se.fit]

Produce Plot of GAM and mean LMER line

points_wavylines_bray_total_year_reduced_gam_nolmer <- ggplot() +
 # geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.2) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r,cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 color = Survey_Name_Season), alpha = 0.5, size = 1) +
    geom_line(data = na.omit(year_survey_unit_expand.dt,cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season)) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt,cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
 # geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0,0.5)) +
  xlab("Year") +
ylab("Total jaccard dissimilarity") +
  theme(legend.position = "null")

points_wavylines_bray_total_year_reduced_gam_nolmer

ggsave(points_wavylines_bray_total_year_reduced_gam_nolmer, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_nolmer.jpg", height = 5, width = 5, unit = "in")


#with lmer

points_wavylines_bray_total_year_reduced_gam <- ggplot() +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0,0.5)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam

ggsave(points_wavylines_bray_total_year_reduced_gam, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam.jpg", height = 6, width = 6, unit = "in")


#ALT
#plot all, but same color scheme (grey)
points_wavylines_bray_total_year_reduced_gam_greyscale <- ggplot() +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = distances_dissimilarities_allyears.r,
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = year_survey_unit_expand.dt,
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = year_survey_unit_expand.dt, aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  rep("black", times = length(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))), name = "Survey Unit") +
  scale_fill_manual(values =  rep("black", times = length(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))), guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year])),
       y = c(0,0.5)) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_greyscale

ggsave(points_wavylines_bray_total_year_reduced_gam_greyscale, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_greyscale.jpg", height = 6, width = 6, unit = "in")

Alternative, color by trend


points_wavylines_bray_total_year_reduced_gam_colorbytrend <- ggplot() +
  geom_ribbon(data = lmer_bray_total_predictions, aes(x = year, ymin = bray_curtis_lmer_preds_lowerCI, ymax = bray_curtis_lmer_preds_upperCI), fill = "grey", alpha = 0.3) +
  geom_point(data = na.omit(distances_dissimilarities_allyears.r, cols = "year_adj"),
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean,
                 fill = Survey_Name_Season), alpha = 0.4, size = 1, shape = 21, color = "white") +
    geom_line(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"),
             aes(x = year,
                 y = bray_glm_mod_fit,
                 color = Survey_Name_Season), alpha = 0.6) +
  geom_ribbon(data = na.omit(year_survey_unit_expand.dt, cols = "year_adj"), aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE, fill =  Survey_Name_Season), alpha=0.1) + #add standard error
  geom_line(data = lmer_bray_total_predictions, aes(x = year, y = bray_curtis_lmer_preds), color = "black") +
    scale_color_manual(values =  color_alpha_order_bytrend, name = "Survey Unit") +
  scale_fill_manual(values =  color_alpha_order_bytrend, guide = "none") +
  theme_classic() +
  lims(x = c(min(distances_dissimilarities_allyears.r[,year]),max(distances_dissimilarities_allyears.r[,year]))) +
  xlab("Year") +
ylab("β-diversity") +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

points_wavylines_bray_total_year_reduced_gam_colorbytrend

ggsave(points_wavylines_bray_total_year_reduced_gam_colorbytrend, path = here::here("figures"), filename ="points_wavylines_bray_total_year_reduced_gam_colorbytrend.jpg", height = 6, width = 6, unit = "in")

#plot each independently for supplement
#all survey names = 
all_survey_names <- sort(unique(distances_dissimilarities_allyears.r$Survey_Name_Season))
#list of plots
points_wavylines_bray_total_year_reduced_gam_individual <- list()
for (i in 1:length(all_survey_names)) {
points_wavylines_bray_total_year_reduced_gam_individual[[i]] <- ggplot() +
  geom_point(data = distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i]],
             aes(x = year,
                 y = bray_curtis_dissimilarity_total_mean), alpha = 0.4, color = "black") +
    geom_line(data = year_survey_unit_expand.dt[Survey_Name_Season == all_survey_names[i]],
             aes(x = year,
                 y = bray_glm_mod_fit), alpha = 0.6) +
  geom_ribbon(data = year_survey_unit_expand.dt[Survey_Name_Season == all_survey_names[i]], aes(x = year, ymin=bray_glm_mod_fit-bray_glm_mod_fit_SE, ymax=bray_glm_mod_fit+bray_glm_mod_fit_SE), alpha=0.1) + #add standard error
  theme_classic() +
#  lims(x = c(min(distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i],year]),max(distances_dissimilarities_allyears.r[Survey_Name_Season == all_survey_names[i],year])),
#       y = c(0,0.5)) +
  xlab("Year") +
ylab("beta-diversity") +
  facet_wrap(~Survey_Name_Season, ncol = 5) +
  theme(legend.position = "null", axis.text = element_text(size = 15), axis.title = element_text(size = 15))

print(points_wavylines_bray_total_year_reduced_gam_individual[[i]])

}

saveRDS(points_wavylines_bray_total_year_reduced_gam_individual, here::here("figures","points_wavylines_bray_total_year_reduced_gam_individual.Rds"))

#print to pdf
library(gridExtra)

ggsave(
   filename = here::here("figures","points_wavylines_bray_total_year_reduced_gam_individual.pdf"), 
   plot = marrangeGrob(points_wavylines_bray_total_year_reduced_gam_individual, nrow=1, ncol=1), 
   width = 8.5, height = 11
)
Warning: Removed 1 rows containing missing values (`geom_point()`).Warning: Removed 1 row containing missing values (`geom_line()`).Warning: no non-missing arguments to max; returning -InfWarning: Removed 1 rows containing missing values (`geom_point()`).Warning: Removed 1 row containing missing values (`geom_line()`).Warning: no non-missing arguments to max; returning -InfWarning: Removed 1 rows containing missing values (`geom_point()`).Warning: Removed 1 row containing missing values (`geom_line()`).Warning: no non-missing arguments to max; returning -Inf

Merge jaccard versus Year plot with GAMS and Region vs. coefficient plot for LMERs

LS0tCnRpdGxlOiAiWWVhciBUb3RhbCBKYWNjYXJkIERpc3NpbWlsYXJpdHkgTW9kZWxzICIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKVGhpcyBjb2RlIGlzIFNjcmlwdCA0ZCBmb3IgS2l0Y2hlbCBldCBhbC4gIkJpb3RpYyBob21vZ2VuaXphdGlvbiwgdGhlIGV4Y2VwdGlvbiBhbmQgbm90IHRoZSBydWxlIGZvciBtYXJpbmUgZmlzaCBjb21tdW5pdGllcyIgbWFudXNjcmlwdC4KCi0gVGhpcyBwcm9qZWN0IGlzIGEgY29sbGFib3JhdGl2ZSBlZmZvcnQgdG8gZGVzY3JpYmUgY2hhbmdlcyBpbiB0YXhvbm9taWMgY29tcG9zaXRpb24gIG9mIGZpc2ggY29tbXVuaXRpZXMgYXJvdW5kIHRoZSB3b3JsZC0tYXMgc2FtcGxlZCBieSBib3R0b20gdHJhd2wgc3VydmV5cy4KCi0gQ29kZSBieSBab8OrIEouIEtpdGNoZWwKClNFU1NJT04gSU5GTyBUTyBETwoKYGBge3Igc2V0dXB9CmxpYnJhcnkoZGF0YS50YWJsZSkKbGlicmFyeSh2ZWdhbikKbGlicmFyeShzZikKbGlicmFyeShjb25jYXZlbWFuKSAjcG9seWdvbiBhcm91bmQgcG9pbnRzCmxpYnJhcnkoYmV0YXBhcnQpICNhbGxvd3MgdXMgdG8gcGFydGl0aW9uIGJldGEgZGl2ZXJzaXR5CmxpYnJhcnkoZ2Vvc3BoZXJlKQpsaWJyYXJ5KGdncHVicikgI3N0YXRfcmVnbGluZV9lcXVhdGlvbgpsaWJyYXJ5KG5sbWUpCmxpYnJhcnkobWdjdikgI3RvIG1ha2UgZ2FtCmxpYnJhcnkoY293cGxvdCkKbGlicmFyeShsbWU0KQoKI1B1bGwgRGlzc2ltaWxhcml0eSBNZWFucwpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIgPC0gcmVhZFJEUyhoZXJlOjpoZXJlKCJvdXRwdXQiLCJkaXNzaW1pbGFyaXRpZXMiLCAiZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yLnJkcyIpKQoKI21ha2Ugc3VydmV5IGFuZCBzdXJ2ZXkgdW5pdCBmYWN0b3JzCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssc3VydmV5Oj1mYWN0b3Ioc3VydmV5KV1bLHN1cnZleV91bml0Oj1mYWN0b3Ioc3VydmV5X3VuaXQpXQoKI2FkanVzdCB5ZWFycwpkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJfYWRqIDo9IHllYXItbWluKHllYXIpKzFdCgojYWRkIG5ldyB2YXJpYWJsZSBmb3IgeWVhciBpbiBzZXF1ZW5jZSBwZXIgcmVnaW9uCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssZmlyc3RfeWVhciA6PSBtaW4oeWVhciksLihzdXJ2ZXlfdW5pdCldCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssbGFzdF95ZWFyIDo9IG1heCh5ZWFyKSwuKHN1cnZleV91bml0KV0KCiNkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJfaW5fc2VxIDo9IHllYXItZmlyc3RfeWVhcisxXQoKZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyc19zYW1wbGVkIDo9IGxhc3RfeWVhci1maXJzdF95ZWFyKzFdCgoKYGBgCgojIyNQYWxldHRlIGZvciBQbG90dGluZwpQYWxldHRlIGZvciBwbG90dGluZyBhbGwgMzcgc3VydmV5IHVuaXRzCmBgYHtyIGxpbmsgY29sb3JzIHRvIHN1cnZleSB1bml0c30Kc3VydmV5X3VuaXQubGlzdCA8LSBsZXZlbHMoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxzdXJ2ZXlfdW5pdF0pCgpwYWxldHRlXzQyIDwtIGMoCiAgIiM1QTUxNTYiLCAjQUkKICAiI0RGMDBEQiIsICNCSVRTLTEKICAiI0RCOEVEQSIsICNCSVRTLTQKICAiI0Y2MjIyRSIsICNDSEwKICAiI0Y4QTE5RiIsICNERk8tTkYKICAiIzE2RkYzMiIsICNERk8tUUNTCiAgIiMzMjVBOUIiLCAjRUJTCiAgIiMzMjgzRkUiLCAjRVZIT0UKICAiI0ZFQUYxNiIsICNGUi1DR0ZTCiAgIiNmY2NiNmQiLCAjR01FWC1GYWxsCiAgIiMxQzgzNTYiLCAjR01FWC1TdW1tZXIKICAiI0M0NDUxQyIsICNHT0EKICAiIzg1NjYwRCIsICNHUkwtREUKICAiI0IwMDA5RiIsICNHU0wtTgogICIjQkY3OUI4IiwgI0dTTC1TCiAgIiMxQ0JFNEYiLCAjSUNFLUdGUwogICIjNzgyQUI2IiwgI0lFLUlHRlMKICAiIzkwQUQxQyIsICNNRURJVFMKICAiIzZCMDAzQSIsICNOQU0KICAiI0E3NUIwMCIsICNORVVTLUZhbGwKICAiI0UzQjA3MiIsICNORVVTLVNwcmluZwogICIjMDJFOEI2IiwgI05JR0ZTLTEKICAiIzk3RTdENSIsICNOSUdGUy00CiAgIiNCMDAwNjgiLCAjTm9yLUJUUy0zCiAgIiMwMEI5RTMiLCAjTlMtSUJUUy0xCiAgIiM5NUUyRjQiLCAjTlMtSUJUUy0zCiAgIiNCM0NFNzMiLCAjTlotQ0hBVAogICIjNjg5NTAwIiwgI05aLUVDU0kKICAiIzM2NGQwMiIsI05aLVNVQkEKICAiI0FBRjQwMCIsICNOWi1XQ1NJCiAgIiNBQTBERkUiLCAjUFQtSUJUUwogICIjN2Y5ZWI4IiwgI1JPQ0tBTEwKICAiI0ZBMDA4NyIsICNTLUdFT1JHCiAgIiNERUEwRkQiLCAjU0NTLVN1bW1lcgogICIjRkNFRjg4IiwgI1NFVVMtZmFsbAogICIjQTU5NDA1IiwgI1NFVVMtc3ByaW5nCiAgIiNGQ0UxMDAiLCAjU0VVUy1zdW1tZXIKICAiIzU0NDU2MyIsICNTV0MtSUJUUy0xCiAgIiNhMzdmYzciLCAjU1dDLUlCVFMtNAogICIjQzA3NUE2IiwgI1dDQU5OCiAgIiNCRENERkYiLCAjWkFGLUFUTAogICIjMDAzRUZGIiAgI1pBRi1JTkQKKQoKY29sb3JfbGluayA8LSBkYXRhLnRhYmxlKHN1cnZleV91bml0ID0gc3VydmV5X3VuaXQubGlzdCxoZXggPSBwYWxldHRlXzQyKQpgYGAKCkFkZCBuYW1lcyBmb3IgcGxvdHRpbmcKYGBge3IgYWRkIG5hbWVzIGZvciBwbG90dGluZ30KCm5hbWVfaGVscGVyIDwtIGRhdGEudGFibGUoU3VydmV5X05hbWVfU2Vhc29uID0gYygiQWxldXRpYW4gSXNsYW5kcyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCYWx0aWMgU2VhIFExIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJhbHRpYyBTZWEgUTQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQ2hpbGUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTmV3Zm91bmRsYW5kIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlF1ZWVuIENoYXJsb3R0ZSBTb3VuZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFYXN0ZXJuIEJlcmluZyBTZWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQmF5IG9mIEJpc2NheSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFbmdsaXNoIENoYW5uZWwiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiR3VsZiBvZiBNZXhpY28gU3VtbWVyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkd1bGYgb2YgQWxhc2thIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkdyZWVubGFuZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIEd1bGYgb2YgU3QuIExhd3JlbmNlIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlMgR3VsZiBvZiBTdC4gTGF3cmVuY2UiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiSWNlbGFuZCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJJcmlzaCBTZWEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTWVkaXRlcnJhbmVhbiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOYW1pYmlhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5FIFVTIEZhbGwiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTkUgVVMgU3ByaW5nIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk4gSXJlbGFuZCBRMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOIElyZWxhbmQgUTQiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQmFyZW50cyBTZWEgTm9yd2F5IFEzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk4gU2VhIFExIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk4gU2VhIFEzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkNoYXRoYW0gUmlzZSBOWiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJFIENvYXN0IFMgSXNsYW5kIE5aIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlcgQ29hc3QgUyBJc2xhbmQgTloiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUG9ydHVnYWwiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUyBHZW9yZ2lhIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTY290aWFuIFNoZWxmIFN1bW1lciIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0UgVVMgRmFsbCIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0UgVVMgU3ByaW5nIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTRSBVUyBTdW1tZXIiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlcgQ29hc3QgVVMiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkF0bGFudGljIE9jZWFuIFpBIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJJbmRpYW4gT2NlYW4gWkEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJSb2NrYWxsIFBsYXRlYXUiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNjb3RsYW5kIFNoZWxmIFNlYSBRMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU2NvdGxhbmQgU2hlbGYgU2VhIFE0IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJGYWxrbGFuZCBJc2xhbmRzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHdWxmIG9mIE1leGljbyBGYWxsIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCYXJlbnRzIFNlYSBOb3J3YXkgUTEiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlN1Yi1BcmN0aWMgTloiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNjb3RpYW4gU2hlbGYgU3ByaW5nIiksCiAgICAgICAgICAgICAgICAgICAgICAgICAgc3VydmV5X3VuaXQgPSBjKAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkFJIiwgICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJJVFMtMSIsICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJJVFMtNCIsICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkNITCIsICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkRGTy1ORiIsICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkRGTy1RQ1MiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkVCUyIsICAgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkVWSE9FIiwgICAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkZSLUNHRlMiLCAgIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkdNRVgtU3VtbWVyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHT0EiLCAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHUkwtREUiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHU0wtTiIsICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHU0wtUyIsICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJJQ0UtR0ZTIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJJRS1JR0ZTIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJNRURJVFMiLCAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJOQU0iLCAgICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJORVVTLUZhbGwiLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJORVVTLVNwcmluZyIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTklHRlMtMSIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTklHRlMtNCIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTm9yLUJUUy0zIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTlMtSUJUUy0xIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTlMtSUJUUy0zIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTlotQ0hBVCIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTlotRUNTSSIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTlotV0NTSSIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUFQtSUJUUyIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiUy1HRU9SRyIsICAgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0NTLVNVTU1FUiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0VVUy1mYWxsIiwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0VVUy1zcHJpbmciLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIlNFVVMtc3VtbWVyIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJXQ0FOTiIsICAgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJaQUYtQVRMIiwgICAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJaQUYtSU5EIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJST0NLQUxMIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTV0MtSUJUUy0xIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJTV0MtSUJUUy00IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJGQUxLIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJHTUVYLUZhbGwiLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIk5vci1CVFMtMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiTlotU1VCQSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiU0NTLVNQUklORyIKICAgICAgICAgICAgICAgICAgICAgICAgICApKQoKCmNvbG9yX2xpbmsgPC0gY29sb3JfbGlua1tuYW1lX2hlbHBlciwgb24gPSAic3VydmV5X3VuaXQiXQoKYGBgCgoKCiMjTWFrZSBHQU1zCgpKYWNjYXJkCmBgYHtyIGphY2NhcmQgZ2Ftc30KYnJheV9jdXJ0aXNfdG90YWxfZ2FtIDwtIGdhbShicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4gfiB5ZWFyICsgcyhzdXJ2ZXlfdW5pdCwgeWVhciwgYnMgPSAiZnMiLCBtID0gMSksI3JhbmRvbSBzbW9vdGgKICAgICAgICAgICAgICAgICAgICAgICAgICAgIGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIpCgpgYGAKCgojI01ha2UgTE1FUlMKCkJyYXkKKlRoZXNlIGFsbCBjb252ZXJnZWQqCmBgYHtyIGJyYXl9CiNydW5uaW5nIHdpdGggbG1lIGluc3RlYWQgb2YgbG1lciBnYXZlIHNhbWUgcmVzdWx0cywgYnV0IGFsbG93ZWQgZm9yIGNhbGN1bGF0aW9uIG9mIHAtdmFsdWUKYnJheV9jdXJ0aXNfdG90YWxfbG1lIDwtIGxtZShicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4gfiB5ZWFyX2FkaiwgcmFuZG9tID0gKH4xICsgeWVhcl9hZGp8c3VydmV5X3VuaXQpLGRhdGEgPSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIpCgojYnV0IGFsc28gcnVuIHdpdGggbG1lciBmb3IgY29uZmludApicmF5X2N1cnRpc190b3RhbF9sbWVyIDwtIGxtZXIoYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuIH4geWVhcl9hZGogKyAoMSArIHllYXJfYWRqfHN1cnZleV91bml0KSxkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yKQoKc3VtbWFyeShicmF5X2N1cnRpc190b3RhbF9sbWUpCmFub3ZhKGJyYXlfY3VydGlzX3RvdGFsX2xtZSkKCmJyYXlfY3VydGlzX3RvdGFsX2NvZWZzIDwtIGRhdGEudGFibGUoY29lZihicmF5X2N1cnRpc190b3RhbF9sbWUpKQpicmF5X2N1cnRpc190b3RhbF9jb2Vmc1ssc3VydmV5X3VuaXQgOj0gcm93bmFtZXMoY29lZihicmF5X2N1cnRpc190b3RhbF9sbWUpKV1bLFllYXIgOj0gcm91bmQoeWVhcl9hZGosNSldWyxJbnRlcmNlcHQgOj0gcm91bmQoYChJbnRlcmNlcHQpYCwyKV0KVmlldyhicmF5X2N1cnRpc190b3RhbF9jb2VmcykKCmJyYXlfY3VydGlzX3RvdGFsX2NvZWZzIDwtIGJyYXlfY3VydGlzX3RvdGFsX2NvZWZzW2NvbG9yX2xpbmssIG9uID0gInN1cnZleV91bml0Il0KCmJyYXlfY3VydGlzX3RvdGFsX2NvZWZzLmV4cCA8LSBicmF5X2N1cnRpc190b3RhbF9jb2Vmc1ssLihTdXJ2ZXlfTmFtZV9TZWFzb24sIEludGVyY2VwdCwgWWVhcildCgojZXhwb3J0IHRoaXMgdGFibGUKZndyaXRlKGJyYXlfY3VydGlzX3RvdGFsX2NvZWZzLmV4cCwgZmlsZSA9IGhlcmU6OmhlcmUoIm91dHB1dCIsImJyYXlfY3VydGlzX3RvdGFsX2NvZWZzLmV4cC5jc3YiKSkKYGBgCgpHZXQgTE1FUiBtb2RlbCBhcyBwcmVkaWN0aW9ucwpgYGB7cn0KCiMgbmVlZCB0byBzb3J0IG91dCB5ZWFyIGluIHNlcSB2ZXJzdXMgb3ZlcmFsbCB5ZWFyIG1vZGVscwojbmV3IGRhdGEgZm9yIGxtZXIKbG1lcl95ZWFyIDwtIHNlcShtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSksIG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSwgYnkgPSAwLjEpCgpsbWVyX3llYXJfYWRqIDwtIHNlcShtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyX2Fkal0pLCBtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyX2Fkal0pLCBieSA9IDAuMSkKCiNwcmVkaWN0IGF2ZXJhZ2UgbG1lcgpsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMgPC0gZGF0YS50YWJsZSh5ZWFyID0gbG1lcl95ZWFyLCB5ZWFyX2FkaiA9IGxtZXJfeWVhcl9hZGopCgojY29uZmlkZW5jZSBpbnRlcnZhbHMKYnJheV9jdXJ0aXNfdG90YWxfbG1lcl9jb25maW50IDwtIGNvbmZpbnQoYnJheV9jdXJ0aXNfdG90YWxfbG1lcikKCiNwb3B1bGF0ZSBkYXRhIHRhYmxlIG9mIGxtZXIgcHJlZGljdGlvbnMKbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzIDo9IGZpeGVmKGJyYXlfY3VydGlzX3RvdGFsX2xtZXIpW1sxXV0gKyBmaXhlZihicmF5X2N1cnRpc190b3RhbF9sbWVyKVtbMl1dICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX2xvd2VyQ0kgOj0gYnJheV9jdXJ0aXNfdG90YWxfbG1lcl9jb25maW50WzVdICsgYnJheV9jdXJ0aXNfdG90YWxfbG1lcl9jb25maW50WzZdICogeWVhcl9hZGpdWyxicmF5X2N1cnRpc19sbWVyX3ByZWRzX3VwcGVyQ0kgOj0gYnJheV9jdXJ0aXNfdG90YWxfbG1lcl9jb25maW50WzExXSArIGJyYXlfY3VydGlzX3RvdGFsX2xtZXJfY29uZmludFsxMl0gKiB5ZWFyX2Fkal0KYGBgCgoKCiMjI01vdmUgZm9yd2FyZCB3aXRoIEphY2NhcmQgVG90YWwgKG90aGVycyB0byBzdXBwbGVtZW50KQoKCkNvZWZmaWNpZW50cyBmb3IgTE1FUiBieSBzdXJ2ZXlfdW5pdApgYGB7cn0KI3VuaXF1ZSBzdXJ2ZXkgdW5pdCB5ZWFyIGNvbWJvcwpzdXJ2ZXlfdW5pdF9zYW1wbGluZ195ZWFycyA8LSB1bmlxdWUoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWywuKHN1cnZleV91bml0LCB5ZWFyX2FkaiwgeWVhciwgeWVhcnNfc2FtcGxlZCldKQoKIyBzZWUgZ3JvdXAgY29lZmZpY2llbnRzCm1vZGVsX2NvZWZzX3JlZHVjZWQgPC0gZGF0YS50YWJsZSh0cmFuc2Zvcm0oYXMuZGF0YS5mcmFtZShyYW5lZihicmF5X2N1cnRpc190b3RhbF9sbWVyKSksIGx3ciA9IGNvbmR2YWwgLSAxLjk2KmNvbmRzZCwgdXByID0gY29uZHZhbCArIDEuOTYqY29uZHNkKSkKI2h0dHBzOi8vc3RhY2tvdmVyZmxvdy5jb20vcXVlc3Rpb25zLzY5ODA1NTMyL2V4dHJhY3QtdGhlLWNvbmZpZGVuY2UtaW50ZXJ2YWxzLW9mLWxtZXItcmFuZG9tLWVmZmVjdHMtcGxvdHRlZC13aXRoLWRvdHBsb3RyYQoKCiNPTkxZIFNMT1BFUwptb2RlbF9jb2Vmc19yZWR1Y2VkIDwtIG1vZGVsX2NvZWZzX3JlZHVjZWRbdGVybSA9PSAieWVhcl9hZGoiLF0KCm1vZGVsX2NvZWZzX3JlZHVjZWRbLHN1cnZleV91bml0IDo9IGdycF1bLHllYXJfYWRqIDo9IGNvbmR2YWxdCgojbWVyZ2Ugd2l0aCBkdXJhdGlvbiBvZiBzdXJ2ZXkKbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGggPC0gbW9kZWxfY29lZnNfcmVkdWNlZFtzdXJ2ZXlfdW5pdF9zYW1wbGluZ195ZWFycywgb24gPSAic3VydmV5X3VuaXQiXQoKCm1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoWyx5ZWFyc19zYW1wbGVkIDo9IGFzLm51bWVyaWMoeWVhcnNfc2FtcGxlZCldWyxEaXJlY3Rpb25hbF9DaGFuZ2UgOj0gaWZlbHNlKHllYXJfYWRqID4gMCwgIkRpZmZlcmVudGlhdGlvbiIsIkhvbW9nZW5pemF0aW9uIildCgojZG9lcyBpdCBjcm9zcyB6ZXJvPwptb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aFssc2lnbmlmaWNhbnQgOj0gaWZlbHNlKGx3ciA+MCAmIHVwcj4wLFQsaWZlbHNlKGx3cjwwICYgdXByPDAsVCxGKSldCgojc2lnbmlmaWNhbnQgZGlyZWN0aW9uYWwgY2hhbmdlCm1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoWyxEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnIDo9IGlmZWxzZSh5ZWFyX2FkaiA+IDAgJiBzaWduaWZpY2FudCA9PSBULCAiRGlmZmVyZW50aWF0aW9uIixpZmVsc2UoeWVhcl9hZGogPCAwICYgc2lnbmlmaWNhbnQgPT0gVCwgIkhvbW9nZW5pemF0aW9uIiwgIk5vIHRyZW5kIGluXG5kaXNzaW1pbGFyaXR5IikpXQoKCiNtaW4gbWF4IGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXMKbWluX2JyYXlfcmVkdWNlZCA8LSBtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyxicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW5dLCBuYS5ybSA9IFQpCm1heF9icmF5X3JlZHVjZWQgPC0gbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclssYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuXSwgbmEucm0gPSBUKQoKbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGggPC0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGhbY29sb3JfbGluaywgb24gPSAic3VydmV5X3VuaXQiXQoKI2RlbGV0ZSBhbnkgTkFzCm1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoIDwtIG5hLm9taXQobW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGNvbHMgPSAic2lnbmlmaWNhbnQiKQoKI29yZGVyIHRhYmxlIGJ5IGNvZWZmaWNpZW50CnNldG9yZGVyKG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCB5ZWFyX2FkaikKCmphY2NhcmRfdG90YWxfbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgudW5pcXVlIDwtIHVuaXF1ZShtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aFssLihjb25kdmFsLGNvbmRzZCwgbHdyLCB1cHIsIHN1cnZleV91bml0LCB5ZWFyX2FkaiwgeWVhcnNfc2FtcGxlZCwgRGlyZWN0aW9uYWxfQ2hhbmdlLCBoZXgsIFN1cnZleV9OYW1lX1NlYXNvbiwgc2lnbmlmaWNhbnQsIERpcmVjdGlvbmFsX0NoYW5nZV9zaWcpXSkgCgojZXh0cmFjdCBjb2xvciBoZXhlcwojeWVhciBhZGogY29lZiBvcmRlcgpjb2xvcl95ZWFyX2Fkal9vcmRlciA8LSBqYWNjYXJkX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZVssaGV4XQoKI2FscGhhYmV0aWNhbCBvcmRlcgpqYWNjYXJkX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZS5hbHBoYSA8LSBzZXRvcmRlcihqYWNjYXJkX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZSwgU3VydmV5X05hbWVfU2Vhc29uKQoKamFjY2FyZF90b3RhbF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUuYWxwaGFbLHRyZW5kX2NvbG9yIDo9IGlmZWxzZShEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnID09ICJIb21vZ2VuaXphdGlvbiIsICIjZTdhYzViIiwgaWZlbHNlKERpcmVjdGlvbmFsX0NoYW5nZV9zaWcgPT0gIkRpZmZlcmVudGlhdGlvbiIsIiM5MWM4NzQiLCIjY2JiZmRlIikpXQoKY29sb3JfYWxwaGFfb3JkZXIgPC0gamFjY2FyZF90b3RhbF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUuYWxwaGFbLGhleF0KY29sb3JfYWxwaGFfb3JkZXJfYnl0cmVuZCA8LSBqYWNjYXJkX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZS5hbHBoYVssIHRyZW5kX2NvbG9yXQoKc2F2ZVJEUyhqYWNjYXJkX3RvdGFsX21vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLnVuaXF1ZSwgaGVyZTo6aGVyZSgib3V0cHV0IiwicmVnaW9uX3N0YXRzIiwiamFjY2FyZF90b3RhbF9tb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aC51bmlxdWUuUmRzIikpCmBgYAoKQmFyIFBsb3QgQ29lZmZpY2llbnQgTE1FUgpgYGB7ciBiYXIgcGxvdCBvZiBjb2VmZmljaWVudHN9CmphY2NhcmRfVE9UQUxfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWQgPC0gZ2dwbG90KCkgKwogICAgZ2VvbV9lcnJvcmJhcihkYXRhID0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGFlcyh4ID0gcmVvcmRlcihTdXJ2ZXlfTmFtZV9TZWFzb24sIHllYXJfYWRqKSAsIHkgPSB5ZWFyX2FkaiwgbGFiZWwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24sIHltaW4gPSBsd3IsIHltYXggPSB1cHIpLCBmaWxsID0gImdyZXkiLCB3aWR0aCA9IDApICsgI2FkZCBjb25maWRlbmNlIGludGVydmFscwogIGdlb21fcG9pbnQoZGF0YSA9IG1vZGVsX2NvZWZzX3JlZHVjZWRfbGVuZ3RoLCBhZXMoeCA9IHJlb3JkZXIoU3VydmV5X05hbWVfU2Vhc29uLCB5ZWFyX2FkaikgLCB5ID0geWVhcl9hZGosIGxhYmVsID0gU3VydmV5X05hbWVfU2Vhc29uLCBzaXplID0geWVhcnNfc2FtcGxlZCwgZmlsbCA9IERpcmVjdGlvbmFsX0NoYW5nZV9zaWcsIGNvbG9yID0gRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZyksIHN0YXQgPSAnaWRlbnRpdHknLCBzaGFwZSA9IDIxKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYygid2hpdGUiLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIiwgZ3VpZGU9Im5vbmUiKSArCiAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IGMoImJsYWNrIiwiYmxhY2siLCJncmV5IiksIG5hbWUgPSAiRGlzc2ltaWxhcml0eSB0cmVuZCIsIGd1aWRlPSJub25lIikgKwogIHNjYWxlX3NpemVfYmlubmVkKHJhbmdlID0gYygxLDgpLCBuYW1lID0gIlN1cnZleSBwZXJpb2QgbGVuZ3RoIikgKwogIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDApICsKICBzY2FsZV95X2NvbnRpbnVvdXMoYnJlYWtzID0gc2VxKC0wLjAwNSwgMC4wMDc1LCBieSA9IDAuMDAyNSksIGxhYmVscyA9IGMoIi0wLjAwNSIsIiIsIjAiLCAiIiwgIjAuMDA1IiwgICIiKSkgKwogIHhsYWIoIlN1cnZleSB1bml0IikgKwogIHlsYWIoIs6yLWRpdmVyc2l0eSB0cmVuZCIpICsgI3RvdGFsIGphY2NhcmQgZGlzc2ltaWxhcml0eSB0cmVuZAogIGNvb3JkX2ZsaXAoKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICB0aGVtZShheGlzLnRleHQueSA9IGVsZW1lbnRfdGV4dChjb2xvdXIgPSBjb2xvcl95ZWFyX2Fkal9vcmRlciwgZmFjZSA9ICJib2xkIiksIGF4aXMudGl0bGUueSA9IGVsZW1lbnRfYmxhbmsoKSwgYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGxlZ2VuZC5wb3NpdGlvbiA9IGMoMC4yLDAuOCksIGxlZ2VuZC5kaXJlY3Rpb24gPSAidmVydGljYWwiLCBsZWdlbmQuYmFja2dyb3VuZCA9IGVsZW1lbnRfcmVjdChmaWxsID0gInRyYW5zcGFyZW50IikpCgojcHVsbCBsZWdlbmQgZm9yIGhvbW9nZW5pemF0aW9uCmRpcmVjdGlvbmFsX2NoYW5nZV9sZWdlbmRfcGxvdCA8LSBqYWNjYXJkX1RPVEFMX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkICsgCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gYygid2hpdGUiLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIikgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCJibGFjayIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGgiLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJyaWdodCIsIGxlZ2VuZC5iYWNrZ3JvdW5kID0gZWxlbWVudF9yZWN0KGZpbGw9ICJ0cmFuc3BhcmVudCIpKSArCiAgZ3VpZGVzKGNvbG91ciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KHNpemU9NikpKQoKCmphY2NhcmRfVE9UQUxfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWQKCiNBTFQgZ3JleSBzY2FsZQpqYWNjYXJkX1RPVEFMX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkX2dyZXlzY2FsZSA8LSBnZ3Bsb3QoKSArCiAgICBnZW9tX2Vycm9yYmFyKGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgeW1pbiA9IGx3ciwgeW1heCA9IHVwciksIGZpbGwgPSAiZ3JleSIsIHdpZHRoID0gMCkgKyAjYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzCiAgZ2VvbV9wb2ludChkYXRhID0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGFlcyh4ID0gcmVvcmRlcihTdXJ2ZXlfTmFtZV9TZWFzb24sIHllYXJfYWRqKSAsIHkgPSB5ZWFyX2FkaiwgbGFiZWwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24sIHNpemUgPSB5ZWFyc19zYW1wbGVkLCBmaWxsID0gRGlyZWN0aW9uYWxfQ2hhbmdlX3NpZywgY29sb3IgPSBEaXJlY3Rpb25hbF9DaGFuZ2Vfc2lnKSwgc3RhdCA9ICdpZGVudGl0eScsIHNoYXBlID0gMjEpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSBjKCJ3aGl0ZSIsImJsYWNrIiwiZ3JleSIpLCBuYW1lID0gIkRpc3NpbWlsYXJpdHkgdHJlbmQiLCBndWlkZT0ibm9uZSIpICsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiYmxhY2siLCJibGFjayIsImdyZXkiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIiwgZ3VpZGU9Im5vbmUiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGgiKSArCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCkgKwogIHNjYWxlX3lfY29udGludW91cyhicmVha3MgPSBzZXEoLTAuMDA1LCAwLjAwNzUsIGJ5ID0gMC4wMDI1KSwgbGFiZWxzID0gYygiLTAuMDA1IiwiIiwiMCIsICIiLCAiMC4wMDUiLCAgIiIpKSArCiAgeGxhYigiU3VydmV5IHVuaXQiKSArCiAgeWxhYigizrItZGl2ZXJzaXR5IHRyZW5kIikgKyAjdG90YWwgamFjY2FyZCBkaXNzaW1pbGFyaXR5IHRyZW5kCiAgY29vcmRfZmxpcCgpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIHRoZW1lKGF4aXMudGV4dC55ID0gZWxlbWVudF90ZXh0KGZhY2UgPSAiYm9sZCIpLCBheGlzLnRpdGxlLnkgPSBlbGVtZW50X2JsYW5rKCksIGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUueCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBsZWdlbmQucG9zaXRpb24gPSBjKDAuMjUsMC44KSwgbGVnZW5kLmRpcmVjdGlvbiA9ICJ2ZXJ0aWNhbCIpCmBgYApBbHRlcm5hdGl2ZWx5LCB3ZSBjb2xvciAgdGhpcyBwbG90IGJ5IHRyZW5kIGV4cGVyaWVuY2VkCgpgYGB7cn0KIyIjNzNCQTREIiwiI0UwOTYyQyIsIiNjYmJmZGUiCgpqYWNjYXJkX3RvdGFsX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkX2NvbG9yYnl0cmVuZCA8LSBnZ3Bsb3QoKSArCiAgICBnZW9tX2Vycm9yYmFyKGRhdGEgPSBtb2RlbF9jb2Vmc19yZWR1Y2VkX2xlbmd0aCwgYWVzKHggPSByZW9yZGVyKFN1cnZleV9OYW1lX1NlYXNvbiwgeWVhcl9hZGopICwgeSA9IHllYXJfYWRqLCBsYWJlbCA9IFN1cnZleV9OYW1lX1NlYXNvbiwgeW1pbiA9IGx3ciwgeW1heCA9IHVwciksIGZpbGwgPSAiZ3JleSIsIHdpZHRoID0gMCkgKyAjYWRkIGNvbmZpZGVuY2UgaW50ZXJ2YWxzCiAgZ2VvbV9wb2ludChkYXRhID0gbW9kZWxfY29lZnNfcmVkdWNlZF9sZW5ndGgsIGFlcyh4ID0gcmVvcmRlcihTdXJ2ZXlfTmFtZV9TZWFzb24sIHllYXJfYWRqKSAsIHkgPSB5ZWFyX2FkaiwgbGFiZWwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24sIHNpemUgPSB5ZWFyc19zYW1wbGVkLCBjb2xvciA9IERpcmVjdGlvbmFsX0NoYW5nZV9zaWcpLCBzdGF0ID0gJ2lkZW50aXR5JykgKwogIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBjKCIjNzNCQTREIiwiI0UwOTYyQyIsIiNjYmJmZGUiKSwgbmFtZSA9ICJEaXNzaW1pbGFyaXR5IHRyZW5kIiwgZ3VpZGU9Im5vbmUiKSArCiAgc2NhbGVfc2l6ZV9iaW5uZWQocmFuZ2UgPSBjKDEsOCksIG5hbWUgPSAiU3VydmV5IHBlcmlvZCBsZW5ndGhcbih5ZWFycykiKSArCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCkgKwogIHNjYWxlX3lfY29udGludW91cyhicmVha3MgPSBzZXEoLTAuMDA1LCAwLjAwNzUsIGJ5ID0gMC4wMDI1KSwgbGFiZWxzID0gYygiLTAuMDA1IiwiIiwiMCIsICIiLCAiMC4wMDUiLCAgIiIpKSArCiAgeWxpbSgtMC4wMTUsMC4wMDYpICsKICB4bGFiKCJTdXJ2ZXkgdW5pdCIpICsKICB5bGFiKCLOsi1kaXZlcnNpdHkgdHJlbmQiKSArICN0b3RhbCBqYWNjYXJkIGRpc3NpbWlsYXJpdHkgdHJlbmQKICBjb29yZF9mbGlwKCkgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgdGhlbWUoYXhpcy50ZXh0LnkgPSBlbGVtZW50X3RleHQoZmFjZSA9ICJib2xkIiksIGF4aXMudGl0bGUueSA9IGVsZW1lbnRfYmxhbmsoKSwgYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZS54ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGxlZ2VuZC5wb3NpdGlvbiA9IGMoMC4yNSwwLjgpLCBsZWdlbmQuZGlyZWN0aW9uID0gInZlcnRpY2FsIiwgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE0KSwgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGxlZ2VuZC5iYWNrZ3JvdW5kID0gZWxlbWVudF9yZWN0KGZpbGwgPSAidHJhbnNwYXJlbnQiKSkKCmRpcmVjdGlvbmFsX2NoYW5nZV9sZWdlbmRfcGxvdF9jb2xvcmJ5dHJlbmQgPC0gamFjY2FyZF90b3RhbF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZF9jb2xvcmJ5dHJlbmQgKyAKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAicmlnaHQiLCBsZWdlbmQuYmFja2dyb3VuZCA9IGVsZW1lbnRfcmVjdChmaWxsPSAidHJhbnNwYXJlbnQiKSwgCiAgICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNCksIGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKSArCiAgZ3VpZGVzKGNvbG91ciA9IGd1aWRlX2xlZ2VuZChvdmVycmlkZS5hZXMgPSBsaXN0KHNpemU9NikpLCBzaXplID0gIm5vbmUiKQpgYGAKCldhdnkgTGluZSBQbG90IGZvciBHQU1zCgpHZW5lcmF0ZSBwcmVkaWN0ZWQgdmFsdWVzCmBgYHtyIGdlbmVyYXRlIHByZWRpY3RlZCB2YWx1ZXMgR0FNfQoKI2FkZCBjb2xvcnMgYW5kIG5hbWVzIHRvIGZ1bGwgZGlzc2ltaWxhcml0eSBkYXRhIHRhYmxlCmRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciA8LSBkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbY29sb3JfbGluaywgb24gPSAic3VydmV5X3VuaXQiXQoKI2dlbmVyYXRlIG5ldyBkYXRhIHRvIHNtb290aCBsaW5lcyAobmVlZCB5ZWFyIGFuZCBzZWFzb24gc3VydmV5IGNvbWJpbmF0aW9ucykKeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQgPC0gZGF0YS50YWJsZShzdXJ2ZXlfdW5pdCA9IGFzLmNoYXJhY3RlcihOVUxMKSwgeWVhciA9IGFzLm51bWVyaWMoTlVMTCksIHllYXJfYWRqID0gYXMubnVtZXJpYyhOVUxMICkpCgpmb3IgKGkgaW4gMTpsZW5ndGgoc3VydmV5X3VuaXQubGlzdCkpIHsKICAjZ2VuZXJhdGUgeWVhciB2ZWN0b3JzCiAgc3VydmV5X3VuaXRfeWVhcnMgPC0gdW5pcXVlKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucltzdXJ2ZXlfdW5pdCA9PSBzdXJ2ZXlfdW5pdC5saXN0W2ldLC4oc3VydmV5X3VuaXQsIHllYXIsIHllYXJfYWRqKV0pCiAgCiAgeWVhcnMgPC0gc2VxKG1pbihzdXJ2ZXlfdW5pdF95ZWFyc1sseWVhcl0pLCBtYXgoc3VydmV5X3VuaXRfeWVhcnNbLHllYXJdKSwgYnkgPSAwLjEpCiAgCiAgeWVhcl9hZGp1c3QgPC0gc2VxKG1pbihzdXJ2ZXlfdW5pdF95ZWFyc1sseWVhcl9hZGpdKSwgbWF4KHN1cnZleV91bml0X3llYXJzWyx5ZWFyX2Fkal0pLCBieSA9IDAuMSkKICAKICB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdF9hZGRpdGlvbiA8LSBkYXRhLnRhYmxlKHN1cnZleV91bml0ID0gc3VydmV5X3VuaXQubGlzdFtpXSwgeWVhciA9IHllYXJzLCB5ZWFyX2FkaiA9IHllYXJfYWRqdXN0KQogIAogIHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0IDwtIHJiaW5kKHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdF9hZGRpdGlvbikKfQoKI2FkZCBjb2xvcnMgYW5kIG5hbWVzIHRvIGZ1bGwgeWVhciBhbmQgc3VydmV5IHVuaXQgY29tYmluYXRpb24gdGFibGUKeWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQgPC0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHRbY29sb3JfbGluaywgb24gPSAic3VydmV5X3VuaXQiXQpgYGAKCgoKR2V0IG1vZGVsIGFzIHByZWRpY3Rpb25zCmBgYHtyfQojZm9yIHBsb3R0aW5nLCBnZXQgbW9kZWwgYXMgcHJlZGljdGlvbnMKYnJheV9jdXJ0aXNfdG90YWxfZ2FtX3ByZWRpY3Rpb25zIDwtIHByZWRpY3QoYnJheV9jdXJ0aXNfdG90YWxfZ2FtLCBzZS5maXQgPSBUUlVFLCBuZXdkYXRhID0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQpCgojbWVyZ2UgaW50byB0YWJsZQp5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFssYnJheV9nbG1fbW9kX2ZpdCA6PSBicmF5X2N1cnRpc190b3RhbF9nYW1fcHJlZGljdGlvbnMkZml0XVssYnJheV9nbG1fbW9kX2ZpdF9TRSA6PSBicmF5X2N1cnRpc190b3RhbF9nYW1fcHJlZGljdGlvbnMkc2UuZml0XQoKYGBgCgoKUHJvZHVjZSBQbG90IG9mIEdBTSBhbmQgbWVhbiBMTUVSIGxpbmUKYGBge3IgcGxvdCBHQU0gYW5kIG1lYW4gTE1FUiBsaW5lc30KcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fbm9sbWVyIDwtIGdncGxvdCgpICsKICMgZ2VvbV9yaWJib24oZGF0YSA9IGxtZXJfYnJheV90b3RhbF9wcmVkaWN0aW9ucywgYWVzKHggPSB5ZWFyLCB5bWluID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc19sb3dlckNJLCB5bWF4ID0gYnJheV9jdXJ0aXNfbG1lcl9wcmVkc191cHBlckNJKSwgZmlsbCA9ICJncmV5IiwgYWxwaGEgPSAwLjIpICsKICBnZW9tX3BvaW50KGRhdGEgPSBuYS5vbWl0KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMucixjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfY3VydGlzX2Rpc3NpbWlsYXJpdHlfdG90YWxfbWVhbiwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC41LCBzaXplID0gMSkgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9nbG1fbW9kX2ZpdCwKICAgICAgICAgICAgICAgICBjb2xvciA9IFN1cnZleV9OYW1lX1NlYXNvbikpICsKICBnZW9tX3JpYmJvbihkYXRhID0gbmEub21pdCh5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdCxjb2xzID0gInllYXJfYWRqIiksIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAjIGdlb21fbGluZShkYXRhID0gbG1lcl9icmF5X3RvdGFsX3ByZWRpY3Rpb25zLCBhZXMoeCA9IHllYXIsIHkgPSBicmF5X2N1cnRpc19sbWVyX3ByZWRzKSwgY29sb3IgPSAiYmxhY2siKSArCiAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyLCBuYW1lID0gIlN1cnZleSBVbml0IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9ICBjb2xvcl9hbHBoYV9vcmRlciwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSksCiAgICAgICB5ID0gYygwLDAuNSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCJUb3RhbCBqYWNjYXJkIGRpc3NpbWlsYXJpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fbm9sbWVyCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fbm9sbWVyLCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX25vbG1lci5qcGciLCBoZWlnaHQgPSA1LCB3aWR0aCA9IDUsIHVuaXQgPSAiaW4iKQoKI3dpdGggbG1lcgoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW0gPC0gZ2dwbG90KCkgKwogIGdlb21fcmliYm9uKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeW1pbiA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfbG93ZXJDSSwgeW1heCA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfdXBwZXJDSSksIGZpbGwgPSAiZ3JleSIsIGFscGhhID0gMC4zKSArCiAgZ2VvbV9wb2ludChkYXRhID0gbmEub21pdChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIsIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLAogICAgICAgICAgICAgICAgIGZpbGwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNCwgc2l6ZSA9IDEsIHNoYXBlID0gMjEsIGNvbG9yID0gIndoaXRlIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNikgKwogIGdlb21fcmliYm9uKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBjb2xzID0gInllYXJfYWRqIiksIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAgZ2VvbV9saW5lKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeSA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHMpLCBjb2xvciA9ICJibGFjayIpICsKICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSAgY29sb3JfYWxwaGFfb3JkZXIsIG5hbWUgPSAiU3VydmV5IFVuaXQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gIGNvbG9yX2FscGhhX29yZGVyLCBndWlkZSA9ICJub25lIikgKwogIHRoZW1lX2NsYXNzaWMoKSArCiAgbGltcyh4ID0gYyhtaW4oZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSksbWF4KGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pKSwKICAgICAgIHkgPSBjKDAsMC41KSkgKwogIHhsYWIoIlllYXIiKSArCnlsYWIoIs6yLWRpdmVyc2l0eSIpICsKICB0aGVtZShsZWdlbmQucG9zaXRpb24gPSAibnVsbCIsIGF4aXMudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpLCBheGlzLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSkpCgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbQoKZ2dzYXZlKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtLCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtLmpwZyIsIGhlaWdodCA9IDYsIHdpZHRoID0gNiwgdW5pdCA9ICJpbiIpCgojQUxUCiNwbG90IGFsbCwgYnV0IHNhbWUgY29sb3Igc2NoZW1lIChncmV5KQpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9ncmV5c2NhbGUgPC0gZ2dwbG90KCkgKwogIGdlb21fcmliYm9uKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeW1pbiA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfbG93ZXJDSSwgeW1heCA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfdXBwZXJDSSksIGZpbGwgPSAiZ3JleSIsIGFscGhhID0gMC4zKSArCiAgZ2VvbV9wb2ludChkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yLAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2N1cnRpc19kaXNzaW1pbGFyaXR5X3RvdGFsX21lYW4sCiAgICAgICAgICAgICAgICAgZmlsbCA9IFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhID0gMC40LCBzaXplID0gMSwgc2hhcGUgPSAyMSwgY29sb3IgPSAid2hpdGUiKSArCiAgICBnZW9tX2xpbmUoZGF0YSA9IHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LAogICAgICAgICAgICAgYWVzKHggPSB5ZWFyLAogICAgICAgICAgICAgICAgIHkgPSBicmF5X2dsbV9tb2RfZml0LAogICAgICAgICAgICAgICAgIGNvbG9yID0gU3VydmV5X05hbWVfU2Vhc29uKSwgYWxwaGEgPSAwLjYpICsKICBnZW9tX3JpYmJvbihkYXRhID0geWVhcl9zdXJ2ZXlfdW5pdF9leHBhbmQuZHQsIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAgZ2VvbV9saW5lKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeSA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHMpLCBjb2xvciA9ICJibGFjayIpICsKICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSAgcmVwKCJibGFjayIsIHRpbWVzID0gbGVuZ3RoKHVuaXF1ZShkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIkU3VydmV5X05hbWVfU2Vhc29uKSkpLCBuYW1lID0gIlN1cnZleSBVbml0IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9ICByZXAoImJsYWNrIiwgdGltZXMgPSBsZW5ndGgodW5pcXVlKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuciRTdXJ2ZXlfTmFtZV9TZWFzb24pKSksIGd1aWRlID0gIm5vbmUiKSArCiAgdGhlbWVfY2xhc3NpYygpICsKICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSxtYXgoZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yWyx5ZWFyXSkpLAogICAgICAgeSA9IGMoMCwwLjUpKSArCiAgeGxhYigiWWVhciIpICsKeWxhYigizrItZGl2ZXJzaXR5IikgKwogIHRoZW1lKGxlZ2VuZC5wb3NpdGlvbiA9ICJudWxsIiwgYXhpcy50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSAxNSksIGF4aXMudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSkKCnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZQoKZ2dzYXZlKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2dyZXlzY2FsZSwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSJwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9ncmV5c2NhbGUuanBnIiwgaGVpZ2h0ID0gNiwgd2lkdGggPSA2LCB1bml0ID0gImluIikKCmBgYAoKQWx0ZXJuYXRpdmUsIGNvbG9yIGJ5IHRyZW5kCmBgYHtyIGNvbG9yIHdhdnkgbGluZXMgYnkgdHJlbmR9Cgpwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9jb2xvcmJ5dHJlbmQgPC0gZ2dwbG90KCkgKwogIGdlb21fcmliYm9uKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeW1pbiA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfbG93ZXJDSSwgeW1heCA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHNfdXBwZXJDSSksIGZpbGwgPSAiZ3JleSIsIGFscGhhID0gMC4zKSArCiAgZ2VvbV9wb2ludChkYXRhID0gbmEub21pdChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIsIGNvbHMgPSAieWVhcl9hZGoiKSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuLAogICAgICAgICAgICAgICAgIGZpbGwgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNCwgc2l6ZSA9IDEsIHNoYXBlID0gMjEsIGNvbG9yID0gIndoaXRlIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBjb2xzID0gInllYXJfYWRqIiksCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQsCiAgICAgICAgICAgICAgICAgY29sb3IgPSBTdXJ2ZXlfTmFtZV9TZWFzb24pLCBhbHBoYSA9IDAuNikgKwogIGdlb21fcmliYm9uKGRhdGEgPSBuYS5vbWl0KHllYXJfc3VydmV5X3VuaXRfZXhwYW5kLmR0LCBjb2xzID0gInllYXJfYWRqIiksIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFLCBmaWxsID0gIFN1cnZleV9OYW1lX1NlYXNvbiksIGFscGhhPTAuMSkgKyAjYWRkIHN0YW5kYXJkIGVycm9yCiAgZ2VvbV9saW5lKGRhdGEgPSBsbWVyX2JyYXlfdG90YWxfcHJlZGljdGlvbnMsIGFlcyh4ID0geWVhciwgeSA9IGJyYXlfY3VydGlzX2xtZXJfcHJlZHMpLCBjb2xvciA9ICJibGFjayIpICsKICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSAgY29sb3JfYWxwaGFfb3JkZXJfYnl0cmVuZCwgbmFtZSA9ICJTdXJ2ZXkgVW5pdCIpICsKICBzY2FsZV9maWxsX21hbnVhbCh2YWx1ZXMgPSAgY29sb3JfYWxwaGFfb3JkZXJfYnl0cmVuZCwgZ3VpZGUgPSAibm9uZSIpICsKICB0aGVtZV9jbGFzc2ljKCkgKwogIGxpbXMoeCA9IGMobWluKGRpc3RhbmNlc19kaXNzaW1pbGFyaXRpZXNfYWxseWVhcnMuclsseWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbLHllYXJdKSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCLOsi1kaXZlcnNpdHkiKSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kCgpnZ3NhdmUocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kLCBwYXRoID0gaGVyZTo6aGVyZSgiZmlndXJlcyIpLCBmaWxlbmFtZSA9InBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2NvbG9yYnl0cmVuZC5qcGciLCBoZWlnaHQgPSA2LCB3aWR0aCA9IDYsIHVuaXQgPSAiaW4iKQpgYGAKCmBgYHtyfQojcGxvdCBlYWNoIGluZGVwZW5kZW50bHkgZm9yIHN1cHBsZW1lbnQKI2FsbCBzdXJ2ZXkgbmFtZXMgPSAKYWxsX3N1cnZleV9uYW1lcyA8LSBzb3J0KHVuaXF1ZShkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnIkU3VydmV5X05hbWVfU2Vhc29uKSkKI2xpc3Qgb2YgcGxvdHMKcG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbCA8LSBsaXN0KCkKZm9yIChpIGluIDE6bGVuZ3RoKGFsbF9zdXJ2ZXlfbmFtZXMpKSB7CnBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2luZGl2aWR1YWxbW2ldXSA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9wb2ludChkYXRhID0gZGlzdGFuY2VzX2Rpc3NpbWlsYXJpdGllc19hbGx5ZWFycy5yW1N1cnZleV9OYW1lX1NlYXNvbiA9PSBhbGxfc3VydmV5X25hbWVzW2ldXSwKICAgICAgICAgICAgIGFlcyh4ID0geWVhciwKICAgICAgICAgICAgICAgICB5ID0gYnJheV9jdXJ0aXNfZGlzc2ltaWxhcml0eV90b3RhbF9tZWFuKSwgYWxwaGEgPSAwLjQsIGNvbG9yID0gImJsYWNrIikgKwogICAgZ2VvbV9saW5lKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtTdXJ2ZXlfTmFtZV9TZWFzb24gPT0gYWxsX3N1cnZleV9uYW1lc1tpXV0sCiAgICAgICAgICAgICBhZXMoeCA9IHllYXIsCiAgICAgICAgICAgICAgICAgeSA9IGJyYXlfZ2xtX21vZF9maXQpLCBhbHBoYSA9IDAuNikgKwogIGdlb21fcmliYm9uKGRhdGEgPSB5ZWFyX3N1cnZleV91bml0X2V4cGFuZC5kdFtTdXJ2ZXlfTmFtZV9TZWFzb24gPT0gYWxsX3N1cnZleV9uYW1lc1tpXV0sIGFlcyh4ID0geWVhciwgeW1pbj1icmF5X2dsbV9tb2RfZml0LWJyYXlfZ2xtX21vZF9maXRfU0UsIHltYXg9YnJheV9nbG1fbW9kX2ZpdCticmF5X2dsbV9tb2RfZml0X1NFKSwgYWxwaGE9MC4xKSArICNhZGQgc3RhbmRhcmQgZXJyb3IKICB0aGVtZV9jbGFzc2ljKCkgKwojICBsaW1zKHggPSBjKG1pbihkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uID09IGFsbF9zdXJ2ZXlfbmFtZXNbaV0seWVhcl0pLG1heChkaXN0YW5jZXNfZGlzc2ltaWxhcml0aWVzX2FsbHllYXJzLnJbU3VydmV5X05hbWVfU2Vhc29uID09IGFsbF9zdXJ2ZXlfbmFtZXNbaV0seWVhcl0pKSwKIyAgICAgICB5ID0gYygwLDAuNSkpICsKICB4bGFiKCJZZWFyIikgKwp5bGFiKCJiZXRhLWRpdmVyc2l0eSIpICsKICBmYWNldF93cmFwKH5TdXJ2ZXlfTmFtZV9TZWFzb24sIG5jb2wgPSA1KSArCiAgdGhlbWUobGVnZW5kLnBvc2l0aW9uID0gIm51bGwiLCBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZSA9IDE1KSwgYXhpcy50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplID0gMTUpKQoKcHJpbnQocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbFtbaV1dKQoKfQpzYXZlUkRTKHBvaW50c193YXZ5bGluZXNfYnJheV90b3RhbF95ZWFyX3JlZHVjZWRfZ2FtX2luZGl2aWR1YWwsIGhlcmU6OmhlcmUoImZpZ3VyZXMiLCJwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9pbmRpdmlkdWFsLlJkcyIpKQoKI3ByaW50IHRvIHBkZgpsaWJyYXJ5KGdyaWRFeHRyYSkKCmdnc2F2ZSgKICAgZmlsZW5hbWUgPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiwicG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbC5wZGYiKSwgCiAgIHBsb3QgPSBtYXJyYW5nZUdyb2IocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1faW5kaXZpZHVhbCwgbnJvdz0xLCBuY29sPTEpLCAKICAgd2lkdGggPSA4LjUsIGhlaWdodCA9IDExCikKCgpgYGAKCk1lcmdlIGphY2NhcmQgdmVyc3VzIFllYXIgcGxvdCB3aXRoIEdBTVMgYW5kIFJlZ2lvbiB2cy4gY29lZmZpY2llbnQgcGxvdCBmb3IgTE1FUnMKCmBgYHtyfQoKamFjY2FyZF9UT1RBTF9HQU1fTE1FUl9tZXJnZV9sZWdlbmQgPC0gZ2dkcmF3KHhsaW0gPSBjKDAsIDQwLjUpLCB5bGltID0gYygwLCAyMSkpICsKICAgIGRyYXdfcGxvdChwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMSwgeSA9IDEsIHdpZHRoID0gMjAsIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoamFjY2FyZF9UT1RBTF9EaXNzaW1pbGFyaXR5X0NvZWZfZXJyb3JiYXJfcmVkdWNlZCArCiAgICAgICAgdGhlbWUobGVnZW5kLmtleS5zaXplID0gdW5pdCgwLjUsICdjbScpLCAjY2hhbmdlIGxlZ2VuZCBrZXkgc2l6ZQogICAgICAgIGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfdGV4dChzaXplPTE2KSwgI2NoYW5nZSBsZWdlbmQgdGl0bGUgZm9udCBzaXplCiAgICAgICAgbGVnZW5kLnRleHQgPSBlbGVtZW50X3RleHQoc2l6ZT0xNCkpLCAjY2hhbmdlIGxlZ2VuZCB0ZXh0IGZvbnQgc2l6ZSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeCA9IDIwLCB5ID0gMSwgd2lkdGggPTE5LCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KGdldF9sZWdlbmQoZGlyZWN0aW9uYWxfY2hhbmdlX2xlZ2VuZF9wbG90ICsgCiAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNSksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTMpKSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMjcsIHkgPSAxMiwgd2lkdGggPSAzLCBoZWlnaHQgPSAyKSArCiAgZ2VvbV90ZXh0KGFlcyh4ID0gMiwgeSA9IDIwLjcpLCBsYWJlbCA9ICgiYS4iKSwgc2l6ZSA9OCwgZm9udGZhY2UgPSAiYm9sZCIpICsKICBnZW9tX3RleHQoYWVzKHggPTIwLCB5ID0gMjAuNyksIGxhYmVsID0gKCJiLiIpLCBzaXplID04LCBmb250ZmFjZSA9ICJib2xkIikKCgpnZ3NhdmUoamFjY2FyZF9UT1RBTF9HQU1fTE1FUl9tZXJnZV9sZWdlbmQsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0gImphY2NhcmRfVE9UQUxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kLnBuZyIsIGhlaWdodCA9IDgsIHdpZHRoID0gMTQsIHVuaXRzID0gImluIikKCiNBTFQgR1JFWSBTQ0FMRQpqYWNjYXJkX1RPVEFMX0dBTV9MTUVSX21lcmdlX2xlZ2VuZF9ncmV5c2NhbGUgPC0gZ2dkcmF3KHhsaW0gPSBjKDAsIDQwLjUpLCB5bGltID0gYygwLCAyMSkpICsKICAgIGRyYXdfcGxvdChwb2ludHNfd2F2eWxpbmVzX2JyYXlfdG90YWxfeWVhcl9yZWR1Y2VkX2dhbV9ncmV5c2NhbGUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeCA9IDEsIHkgPSAxLCB3aWR0aCA9IDIwLCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KGphY2NhcmRfVE9UQUxfRGlzc2ltaWxhcml0eV9Db2VmX2Vycm9yYmFyX3JlZHVjZWRfZ3JleXNjYWxlICsKICAgICAgICB0aGVtZShsZWdlbmQua2V5LnNpemUgPSB1bml0KDAuNSwgJ2NtJyksICNjaGFuZ2UgbGVnZW5kIGtleSBzaXplCiAgICAgICAgbGVnZW5kLnRpdGxlID0gZWxlbWVudF90ZXh0KHNpemU9MTYpLCAjY2hhbmdlIGxlZ2VuZCB0aXRsZSBmb250IHNpemUKICAgICAgICBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTE0KSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB4ID0gMjAsIHkgPSAxLCB3aWR0aCA9IDE5LCBoZWlnaHQgPSAyMCkgKwogICAgZHJhd19wbG90KGdldF9sZWdlbmQoZGlyZWN0aW9uYWxfY2hhbmdlX2xlZ2VuZF9wbG90ICsgCiAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNSksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTMpKSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyNywgeSA9IDEyLCB3aWR0aCA9IDMsIGhlaWdodCA9IDIpICsKICBnZW9tX3RleHQoYWVzKHggPSAyLCB5ID0gMjAuNyksIGxhYmVsID0gKCJhLiIpLCBzaXplID04LCBmb250ZmFjZSA9ICJib2xkIikgKwogIGdlb21fdGV4dChhZXMoeCA9MjAsIHkgPSAyMC43KSwgbGFiZWwgPSAoImIuIiksIHNpemUgPTgsIGZvbnRmYWNlID0gImJvbGQiKQoKZ2dzYXZlKGphY2NhcmRfVE9UQUxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2dyZXlzY2FsZSwgcGF0aCA9IGhlcmU6OmhlcmUoImZpZ3VyZXMiKSwgZmlsZW5hbWUgPSAiamFjY2FyZF9UT1RBTF9HQU1fTE1FUl9tZXJnZV9sZWdlbmRfZ3JleXNjYWxlLnBuZyIsIGhlaWdodCA9IDgsIHdpZHRoID0gMTQsIHVuaXRzID0gImluIikKCiNBTFQgQ09MT1IgQlkgVFJFTkQKamFjY2FyZF90b3RhbF9HQU1fTE1FUl9tZXJnZV9sZWdlbmRfY29sb3JieXRyZW5kIDwtIGdnZHJhdyh4bGltID0gYygwLCA0MC41KSwgeWxpbSA9IGMoMCwgMjEpKSArCiAgICBkcmF3X3Bsb3QocG9pbnRzX3dhdnlsaW5lc19icmF5X3RvdGFsX3llYXJfcmVkdWNlZF9nYW1fY29sb3JieXRyZW5kLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAxLCB5ID0gMSwgd2lkdGggPSAyMCwgaGVpZ2h0ID0gMjApICsKICAgIGRyYXdfcGxvdChqYWNjYXJkX3RvdGFsX0Rpc3NpbWlsYXJpdHlfQ29lZl9lcnJvcmJhcl9yZWR1Y2VkX2NvbG9yYnl0cmVuZCArCiAgICAgICAgdGhlbWUobGVnZW5kLmtleS5zaXplID0gdW5pdCgwLjUsICdjbScpLCAjY2hhbmdlIGxlZ2VuZCBrZXkgc2l6ZQogICAgICAgIyBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNiksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIyBsZWdlbmQudGV4dCA9IGVsZW1lbnRfdGV4dChzaXplPTE0KQogICAgICAgKSwgI2NoYW5nZSBsZWdlbmQgdGV4dCBmb250IHNpemUpLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyMCwgeSA9IDEsIHdpZHRoID0gMTksIGhlaWdodCA9IDIwKSArCiAgICBkcmF3X3Bsb3QoZ2V0X2xlZ2VuZChkaXJlY3Rpb25hbF9jaGFuZ2VfbGVnZW5kX3Bsb3RfY29sb3JieXRyZW5kICsgCiAgICAgIHRoZW1lKGxlZ2VuZC5rZXkuc2l6ZSA9IHVuaXQoMC41LCAnY20nKSwgI2NoYW5nZSBsZWdlbmQga2V5IHNpemUKICAgICAgICBsZWdlbmQudGl0bGUgPSBlbGVtZW50X3RleHQoc2l6ZT0xNiksICNjaGFuZ2UgbGVnZW5kIHRpdGxlIGZvbnQgc2l6ZQogICAgICAgIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemU9MTUpKSksICNjaGFuZ2UgbGVnZW5kIHRleHQgZm9udCBzaXplKQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHggPSAyNi41LCB5ID0gNywgd2lkdGggPSAzLCBoZWlnaHQgPSAyKSArCiAgZ2VvbV90ZXh0KGFlcyh4ID0gMiwgeSA9IDIwLjcpLCBsYWJlbCA9ICgiYS4iKSwgc2l6ZSA9OCwgZm9udGZhY2UgPSAiYm9sZCIpICsKICBnZW9tX3RleHQoYWVzKHggPTIwLCB5ID0gMjAuNyksIGxhYmVsID0gKCJiLiIpLCBzaXplID04LCBmb250ZmFjZSA9ICJib2xkIikKCmdnc2F2ZShqYWNjYXJkX3RvdGFsX0dBTV9MTUVSX21lcmdlX2xlZ2VuZF9jb2xvcmJ5dHJlbmQsIHBhdGggPSBoZXJlOjpoZXJlKCJmaWd1cmVzIiksIGZpbGVuYW1lID0gImphY2NhcmRfdG90YWxfR0FNX0xNRVJfbWVyZ2VfbGVnZW5kX2NvbG9yYnl0cmVuZC5wbmciLCBoZWlnaHQgPSA4LCB3aWR0aCA9IDE0LCB1bml0cyA9ICJpbiIpCgpgYGAKCg==